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Abstract 


We propose a simple preferential attachment model of growing network 
using the complementary probability of Barabasi-Albert (BA) model, i.e., 


n(fci) oc 1 


——. In this network, new nodes are preferentially attached to 


not well connected nodes. Numerical simulations, in perfect agreement with 
the master equation solution, give an exponential degree distribution. This 
suggests that the power law degree distribution is a consequence of prefer¬ 
ential attachment probability together with “rich get richer” phenomena. 
We also calculate the average degree of a target node at time t (< ks{t) >) 
and its fluctuations, to have a better view of the microscopic evolution of 
the network, and we also compare the results with BA model. 


1 Introduction 


In recent years there is a growing interest to study the evolution of complex 
networks and to develop models that reflect certain properties of the real net¬ 
works using some statistical mechanics techniques, graph theory and computer 
simulations m [211311]. One of the most important properties studied in networks 
is the degree distribution of nodes which is the probability P{k) of a node to have 
degree k. We can distinguish three main laws of degree distribution: Poisson law 


k and 7 represents the 


k 

where P{k) = -, power law with P(k) 

kl 

degree exponent, and exponential law with P{k) ~ where c is constant. 

It appears that in nature most networks follow the last two distribution laws 
referred to above. Barabasi-Albert reinvented price’s power law degree distri¬ 
bution network by introducing a simplified model based on both growth and 
preferential attachment. The resulting scale-free network is widely observed in 
variety of systems such as publication citation networks, many social networks, 
protein and gene networks. However, there are other real networks that follow 
an exponential law, for example. Worldwide Marine Transportation Network [5], 
the North American Power Grid Network [3, neural network of the C.elegans 
[7|, and the Email Network at the University of Rovira i Virgili (ENURV) in 
Spain [ 8 ]- 
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The exponential law seems to be the result of growing network by randomly 
adding new nodes and links. On the other hand the power law seems to appear 
when nodes are added to the network one at time and are linked with nodes 
already well connected. 

Many ideas on the formation of networks have been reviewed in the recent years. 
For example, Barabasi, in his earlier work, asserted that the preferential attach¬ 
ment and growth are both required to generate a scale-free network [^. Actually, 
it seems that growth is not necessary for such a purpose |10]. Furthermore, it 
is intuitive that preferential attachment without rich get richer effect does not 
generate a scale-free network m- Krapivsky et al. [12] have studied non linear 
preferential attachement with Ii{ki) oc kj, and they shown that for 7 < 1 the 
mechanism produces a stretched exponential degree distribution. 

Despite many efforts, consistent theory of networks in evolution is still lacking 
and there is not yet a general principle predicting the topology of a formed net¬ 
work. 

Aiming to understand the formation and the evolution of complex networks, 
many models were introduced to investigate the microscopic processes impli¬ 
cated in the resulting network structure. 

In this context we introduce a simple complex network model growing with lin¬ 
ear preferential attachment mechanism and without rich get richer effect.The 
objective is twofold: first to check if the power law degree distribution remain 
in the absence of the rich get richer scenario and, second, to see for eventual 
microscopic differences between scale-free and homogeneous networks. 

2 Degree distribution 

Similarly to the original BA model, our network evolves according to two mecha¬ 
nisms: the growth and the preferential attachment. Nodes entering the network 
prefer to attach to nodes with low degree, then the probability Il{ki) that one 

of the links of a new node connects to node i depends on its degree ki such that 
kj 

Ylj 

In connection with social networks, if we consider the degree of nodes as describ¬ 
ing the wealth of people in a capitalist society, it is known m that we live in 
a world where rich get richer, but what kind of society we will have if there is 
no favors to rich people, and there is instead a continuous subvention to poor 
people?. 

To implement our idea, we start with niQ nodes, each one with m links. At 
every time step we add a new node with m edges that link the new node to m 
different nodes already present in the network. The probability that the new 

ki 

Y.j kj 

t 

malization constant C is deduced from the condition = I, which gives 

i=l 

C = -. t is the time when the last node was created and represents 

t + mo — 1 

also the number of nodes added to the network. 

For this model, the master equation can be written as: 

{t -|- l)P{k, t -|- I) = tP{k, t) + mll{k — I, t)tP{k — I, t) — mli(k, t)tP{k, t) -|- Sk m, 

( 1 ) 


). The nor- 


node is connected to a node i of degree ki is n(A:j) = C{l 


), where C is a normalization constant. 


n(A:,) = C(I 


2 





where 5 is the Kronecker symbol. 

The corresponding stationary equation takes the form: 


{t + l)P{k) = tP{k) + m(l 


k-1 

2mt + mmo 


) 


tP{k - 1) 
t - 1 


— m(l 


k 

2mt + mniQ 


) 


tP{k) 

t-1 

( 2 ) 


+ 


,m: 


where we used kj = 2mt + mrriQ. 






P{k) 


2mt — {k — 1) 
2t + 2mt — k 


P{k-l), 




for k > m, 


2t 

- 2t + 2mt — m ’ 


for k = m. 


The above recurrence relation yields the following solution: 


(3) 


P{k) = < 


2t 


n 


2t + 2mt — m V 2t + 2mt — j 

j=m+l ^ 


2mt — j + 1 


T ) , for fc > m, 


(4) 


2t 


K2t + 2mt — m 


for k = m. 


Although this equation is not in a closed form, numerical estimation of P{k) is 
straightforward as shown in Fig. [TJ 

We also simulate the network with sizes up to re = 2 x 10®, initial number of 
nodes mo = 3 and m = 2. The simulation results strongly support the analytical 
findings (see Fig. [T|). 

We observed in simulations that k remains less than 40 for t = 2.10®, we then 
take t S> j in Eq. m and we obtain 


P{k) 


1 

1 + m 

< 

1 

- 1 + m 


m \k-Tn-l 

1 + mJ 

for k = m. 


for k > m, 


(5) 


After normalization we get the exponential degree distribution P{k) = Ae 

TTi ~\~ 1 

with A = ln(-). The inset in Fig. [1] shows the exponential form of P{k) 

m 

and the excellent agreement between simulations and theoretical results. This 
clearly confirms that the preferential attachment alone is not sufficient to pro¬ 
duce scale-free networks. 


3 Comparison with BA model 

We search for differences between heterogeneous and homogeneous networks by 
comparing our model with the BA model. The degree distribution alone is not 
enough to characterize networks. Computing others microscopic quantities may 
help to have better insight into their evolution and formation. It turns out that 
scale-free network has nodes with important degree (hubs), while random net¬ 
work has no apparent structure. Evaluating the instantaneous average degree 
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Figure 1: Simulation results (circles) for n = 2.10®, m = 2, ttiq = 3, and numer¬ 
ical solution (solid line) of Eq. dH). In the inset we plot the same data in the 
log-linear scale. 


of target node < ks{t) > and its fluctuations, can give quantitative information 
about hubs in the network. In fact, < ks{t) > is somehow related to the instan¬ 
taneous average degree of hubs, because when choosing randomly nodes, hubs 
have more chance to be selected. 

Firstly, we analyze < ks{t) > and < k‘l{t) > in the BA network 

t 

< ks{t) >= ^ U{ki)ki{t) + moU{ko)ko{t), ( 6 ) 

U=l 


where 11 (A:*) = 


ki{t) 


, ti is the time when the node i was created, and 

2mt -|- mmo 
/co(t) is the degree of initial nodes at time t. 

dki{t) 


Solving the mean field equation 


dt 


Inserting the last expression in Eq. 


= mll{ki), we obtain ki{t) = m (— 

V 2ti -I- mo 

we get 


1 

2 


<k.(t)> = + (7) 

ti - 1 

= ™(M2( + mo)+7-o+5pj:^ + 0(i)), (8) 

, , ^ 1 ,11 1 

where 7 is the Euler constant, and a = —I-h ... H-• 

' ’ 2 3 1 + mo 

Good agreement is obtained as shown in Fig.[2)(a) between Eq. ([ 8 ]) and simulation 
results even for the hrst moments of the evolution. < ks{t) > grows indefinitely 
with time and diverges for infinite network (or t ^ 00 ) due to the fact that, in 
heterogeneous networks, hubs are more likely to be selected and linked with new 
nodes. 
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On the other side, the average degree of the network remains finite miis] since 
the majority of nodes have a small degree and the weight of hubs is small. 



Figure 2: (a) Evolution of < ks{t) > in the BA model, the solid line represents 
Eq. Q. (b) Evolution of fluctuations of < ks{t) >, the solid line represents 
Eq. (jlip . Circles joined by dashed lines in both cases are simulations data 
averaged over 20 runs for m = 2, mo = 3. 


The second moment < kg{t) > is written as 


< k‘^{t) > = '^Il{ki)ki{t) +mon(fco)fco(t) 

ti=i 

1*1 1 

« m^(2t + mo)^ ( -3- + mo^ 

{2ti + mo)2 


(9) 

( 10 ) 


Eor large time, f — 


ti=i 


1 \ ^ 
1 \ 2 


= C( 2 ) ~ 2.612, we obtain < kg{t) >« m^\/^(mQ ^ + 


2.612 — b) with b = 1 -3- H-3- + ... H-3-. 

22 32 (1 + mo)2 

Eluctuations of < ks{t) > are given by 


{Aks{t))‘^ =< ks{t)‘^ > — < ks{t) m^ {ttiq ^ + 2.612 — b)V^ — (ln(2t))^ , 

( 11 ) 


which become arbitrary large when time increases sufficiently. 

Simulation data, in accordance with Eq. m (see Fig. [2Kb)), shows the increasing 
tendency of fluctuations in < ks{t) >. This can be explained by the fact 
that the maximum degree in the network k^ax ~ Vt increases [15] faster than 
< ks{t) >~ ln(f) (Eq. ([8|)) and the difference between the two quantities becomes 
greater with time. 
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We now turn to the same analysis in our model. The mean field evolution 
equation for ki{t) gives 

dkjjt) __ m 

dt {2t + mo){t + niQ — 1) t + niQ — 1 

The solution has the form 


ki(t) = m 


t + mo — 1 
2t + mo 


1 

mo-2 


/ tj + mp - 1 

V 2ti + mo 


1 

mQ —2 


A{ti) + A{t) , 


(13) 


, / f'+mp —1 \ "10-2 

where A(t)= [ ^ ^ - dt'. 

t' + mo- I 

The average value of target node < ks{t) > is obtained immediately for any 
time t by substituting Eq. (fTB]) into Eq. ([U]). The resulting equation is solved 
numerically as shown in Fig. [3Ka). 

For large time and taking t 3> mp, we find A{t) ~ 2 ’"o -2 ln(t), ki{t) ~ m^l + 

In —, then 
ti)' 


t 

< A:s(t) > « y ^ ln(t) - In(ti) + (14) 

ti=i 

~ 2m 


This is the value of the node average connectivity of a network growing with 
uniform attachment resulting in an exponential distribution of degree [S]. 

The second moment is obtained by substituting the corresponding expressions 
of Ti{ki) and ki{t) in Eq. ([9|), we get for large time 


< ki{t) >« + 1 ) 


ti=i 


(15) 


t t 

Making the approximations In(tj) ~ tln(f) — t, and E In(tj)^ « tln(t)^ — 

ti=i ti=i 

2tln{t) + 2t — 2, we find < kg{t) >« 5m^. 

Fluctuations are {Aksit))"^ = < kg{t) > — < ks{t) ps m^. This finding, to¬ 
gether with < ks{t) >~ 2m, show that almost all nodes have the same degree 
as illustrated in Fig. [3Kb). The homogeneity of the network can be explained by 
the fact that the preferential attachment used here doesn’t allow the formation 
of hubs, since it neither allows the rich to get richer, nor it enriches the poor. 


4 Conclusion 

In this work, we have introduced a simple model of complex network with a 
preferential attachment criteria and without ’’rich get richer” effect. The net¬ 
work obtained is homogeneous, which demonstrates the crucial role of the ’’rich 
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Figure 3: (a) Evolution of < ks{t) > in our model, the solid line represents 
Eq. ([HI), (b) Evolution of fluctuations of < ks{t) >, the solid line represents the 
numerical solution of Eq. (jl4p and Eq. m- Circles joined by dashed lines in 
both cases are simulations data averaged over 20 runs for m = 2, niQ = 3. 


get richer” in the topology of the network. Giving preferential treatment to the 
least connected nodes is equivalent to use a random attachment probability. In 
terms of social wealth distribution, Pareto principle [13] doesn’t apply and we 
have instead an exponential distribution of income. 

Computing the instantaneous average degree of a target node and its fluctua¬ 
tions provide more information than the usual average degree of the network, in 
particular we show how the average degree of hubs and its fluctuations diverge 
with time in the BA model, and stay finite in our model. 
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